The isotropic XY model on the 
inhomogeneous periodic chain 



J. P. de Lima"'\ T. F. A. Alves^ L. L. Gongalves"' 

Departamento de Fisica Geral, Instituto de Fisica, Universidade de Sao Paulo, 
CP. 66318, 05315-970, Sdo Paulo, SP, Brazil 

^ Departamento de Fisica, Universidade Federal do Ceard, Campus do Pici, 
CP. 6030, 60451-970, Fortaleza, Ceard, Brazil 



Abstract 

The static and dynamic properties of the isotropic XY-model (s = 1/2) on the 
inhomogeneous periodic chain, composed of N segments with n different exchange 
interactions and magnetic moments, in a transverse field h are obtained exactly at 
arbitrary temperatures. The properties are determined by introducing the gener- 
alized Jordan- Wigner transformation and by reducing the problem to a diagonal- 
ization of a finite matrix of n-th order. The diagonalization procedure is discussed 
in detail and the critical behaviour induced by the transverse field, at T = 0, is 
presented. The quantum transitions are determined by analyzing the behaviour of 
the induced magnetization, defined as 0-/n)Y^'^=il^m < S^j m > where is the 
magnetic moment at site m within the segment j, as a function of the field, and the 
critical fields determined exactly. The dynamic correlations, < <S'|^(t)5^, ^, (0) >, 
and the dynamic susceptibility Xq^il^) ^-i^s also obtained at arbitrary temperatures. 
Explicit results are also presented in the limit T = 0, where the critical behaviour 
occurs, for the static susceptibility Xg^(O) ^ ^ function of the transverse field h, and 
for the frequency dependency of dynamic susceptibility Xq^i^)- Also in this limit, 
the transverse time-correlation < SJ^(t)S'J, ^, (0) >, the dynamic and isothermal 

susceptibilities, Xq^i^) and Xt^, are obtained for the transverse field greater or equal 
than the saturation field. 
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1 Introduction 



Models involving inhomogeneous spin chains have been subject of intensive 
study in recent years motivated by various reasons. Amongst those, the most 
remarkable one is the necessity to understand the unusual new properties pre- 
sented by low dimensional magnetic materials [1,2,3,4] at low temperature, 
which are described in terms of its many-body behaviour and its quantum 
transitions [5]. As manifestation of these properties we can mention the ap- 
pearance of magnetization plateaus as functions of the external magnetic field 
[6.7], the existence of an energy gap between the ground state and first excited 
state at zero field [8] and the presence of quantum critical behaviour [4,9,10]. 

The one-dimensional XY model (s = 1/2) introduced by Lieb et al. [11] plays 
an important role in this context since it constitutes one of the few many- 
body problems which can be exactly solved. The most recent results on the 
inhomogeneous anisotropic model have been obtained by Derzhko et a/. [12] 
and are restricted to the thermodynamic properties. A good review of the 
known results is also presented in his work. 

For the isotropic model, the thermodynamic properties have also been ob- 
tained by Derzhko[13] (see references therein for thermodynamic properties 
of the alternating chain), and the study of the dynamics and of the quantum 
critical behaviour has been restricted to the alternating chain [14,15,16], and 
for the alternating superlattice[17] . 

In this paper we will also consider the isotropic XY model in a transverse 
field on the inhomogeneous periodic chain consisting of N unit cells with n 
sites, which has been studied by Derzkho [13]. Within the cells we can have 
n different exchange constants as well as magnetic moments, and the model 
corresponds to an extension of the alternating superlattice. An extensive study 
of the quantum critical behaviour will be presented and the dynamic properties 
in the field direction determined for arbitrary temperature. Within a new 
formalism, we have been able to solve the model exactly, and present new 
and more general results from those presented in the previous papers[17,18]. 
In particular, we have been able to obtain the dynamic correlation in the 
xy-plane, at T = 0, for transverse field greater than the saturation field. 

In section 2 we discuss in detail the diagonalization of the model and present 
the analytical results for the chains composed of cells with two, three and 
four sites. The excitation spectrum is also obtained, by a different method, in 
Appendix A, and it is shown under which conditions the gap in the excitation 
spectrum is suppressed. 
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Fig. 1. Unit cell of the periodic chain. 



The induced magnetization and the isothermal susceptibility Xt obtained 
in section 3 and we also determine the critical fields associated with the quan- 
tum second order phase transitions. The static and dynamic correlations in 
the field direction are obtained in section 4 and, in the section 5, we deter- 
mine the longitudinal dynamic susceptibility Xq^{^)- The dynamic correlation 

(^^1 m{^)^i+i m'{^)) obtained, at T = and for external fields greater than 
the saturation field, in appendix B. Under these conditions, the isothermal and 
dynamic transverse susceptibility x^'^' are obtained in section 6, and finally in 
section 7 we summarize the main results and present the conclusions. 



2 The diagonalization of the Hamiltonian 



We consider the isotropic XY model (s = 1/2) on the inhomogeneous periodic 
chain with cells, n sites per cell, and lattice parameter a, in a transverse 
field, whose unit cell is shown Fig. 1. The Hamiltonian is given by 



n-l 



H — X] I XI l^mhSi j^ + X^ -^TO [Sl,m^l,m+1 + ^l,m^i 
1=1 \m=l m=l 



y 

l,m+l 



+ 



(1) 



where the parameters Ji^m are the exchange coupling between nearest-neighbour, 

the magnetic moments, h the external field and we have assumed periodic 
boundary conditions. If we introduce the ladder operators 



(2) 



and the generalized Jordan- Wigner transformation [19] 
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I l'=lm'=l m'=l ) 

where c;,^ and c],^ are fermion annihilation and creation operators, we can 
write the Hamiltonian as [20] 

H = H+P^ + H-P-, (4) 

where 

= - E 1 E \^rnh ( c[„Q,„ - 9 ) + II ^ (4mCi,m+l + cJ^+iQ,^) [ - 
Z=l lm=l ^ m=l ^ J 

- II Y (4nC/+l,l + cJ+l.lC/.n) ± Y (4,nCl,l + A,\CN,r) , (5) 



and 



- ^. (6) 



with P given by 

(JV n \ 
^^1212 4,mCl,m ■ (7) 
1=1 m=l J 

As it is well known [20,21,22], since the operator P commutes with the Hamil- 
tonian, the eigenstates have definite parity, and P~(P+) corresponds to a 
projector into a state of odd (even) parity. 

Introducing periodic and anti-periodic boundary conditions on c's for H~ and 
respectively, the wave- vectors in the Fourier transform [23], 



Q,m = -j^ H exp {-iqdl) Ag^^, 
1 ^ 

ViV 

are given Q~ — ^ for periodic condition and q'^ — for anti-periodic 

condition, with r = 0, ±1, ....,±N/2, and H~ and H'^ can be written in the 
form 

H^^YHq^^ (9) 

where 
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" / 1 \ 

m=l ^ 



n— 1 T 
m=l ^ 



(10) 



Although i7~ and do not commute, it can be shown that in the thermo- 
dynamic hmit all the static properties of the system can be obtained in terms 
of H~ or . However, even in this limit, some dynamic properties depend 
on H- and H+ [20,21,22]. Since 



[Hq, Hqi] = 0, 



where we make the identification g = we can diagonalize the Hamiltonian 
by introducing the canonical transformation 



^q,m — E '^q,km^q,k, ^q,m ~ E ^q,km^q,k^ 



k=l 



k=l 



and by imposing the condition 

[^q,ki Hq] £q,k^q,ki 

which leads, for the coefficients Uq^km, to the equation 



^Uq^kl^ 

Uq^k2 



\y-q,kn J 



^q,k 



^Uq^kl^ 



\^'^q,kn j 



where is given by 



h 


Jl 

2 





• • • ^ exp ( 


Jl 

2 


h2 


2 








2 


h 


^3 : 

2 






ia 

2 


■•• ■-. 









h 

• ""n— 1 2 



\^ exp {iqd) • • • 



Jn- 



[12) 



(13) 



(14) 



(15) 



and the u's satisfy the orthogonality relations 
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Yl ^q,kmU*q^k'm = ^kk', (16) 
m=l 
n 

Uq.km'U'q.km' = ^mm'- (17) 

fe=l 

Therefore the Hamiltonian can be written in the diagonal form 

H,^E^^A^ikUk-b: (18) 
k 

where the spectrum eg of Hg is determined from the determinantal equation 

det{Aq - 6ql) = 0. (19) 

In passing, we would like to note that for uniform magnetic moments, fim = ^, 
the term — Ydm f^hSi^ commutes with the Hamiltonian, and consequently the 
effect of the field is to shift the spectrum. 

By using (16) we can express the operators ^'s in terms of A's which are given 

by 

n n 

^q,k — X] '^*q,km^q,m : ^g,k = X] '"■9,':"^^g,m ' (^0) 

m=l m=l 

and from eq.(8) we obtain, 

N n 

Uk = ^ exp (iqdl) ulkmCi,m, 

N n 

Cl,k = ^ EE exp i-iqdl) Ug^kmclm, (21) 

ViV 1=1 m=i 

and their inverse, as 

Cl,m = i E E {-iQdl) Uq,kmCq,k, 

ViV q 

4,m = ^ E E (^^^0 K,km^lk- (22) 



^ q k 



The solution of eqs.(14) and (19) can be obtained analytically for n < 4, and 
numerically for n > 4. In particular, the exact dispersion relations for the 
cases n — 2,3,4 are given below. 

For n = 2 is given by 

^q = -^(/^i + l^2)h ± -^(//i - At2)2/i2 + j2 + j2 ^ J2 cos(2g), (23) 
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and for n = 3 by [24] 



£9 = -^(/^i + At2 + Ai3)/i + 2v^cos^^^^^j , (24) 

where p — 0, 2, 4 and 9q — arccos {Rq/ V— -R^) , with R and Rq given by 

3a2 - a? 90201 - 27oq - 2of 
^ = g ' ^1 ' (^5) 



where 



«i = (ah + + /^a) ^, 

02 = (//1//2 + A«lA«3 + A*2A*3) ~ 4 (-^1 + -^2 + Jl) > 

/ii/i2/X3/i^ - ^ (/^i-''2 + /^2-^| + Ai3</i) + ^-^1-^2-^3 cos(3g). (26) 



aq 



For n = 4, the four branches of the dispersion relation are obtained from the 
expression [24] 



A ^JA^+Wq ± J-(A + 2Yq) - 2(^1 T B^/JA^+Wq) 

+ — 7. , (27) 



by considering the sign combinations (+ + +, , — I — , H h), and 

where 



3^2 
'~8" 



-U-2,3 
6 V 3 



cos( ), 



Sn 



arccos 



27_ 



' 12 

3^4 5y42 



^? , AiG, Bl 



108 



256 



+ Dq, 



(28) 



and 
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B = (/X2Atl + A*4AH + /ia/il + /i4A*2 + /^3A«2 + A«4A*3)/i^ - 



4 

,3 



n ,4 (■^1/^3/^2 + ^1/^2/^1 + ^1/^4/^1 + ^1^4/^3)^^ , 
-Dg = Ai4/i3/t2/il/i ^ h 

The excitation spectrum, given by eq.(19), can also be obtained by a transfer 
matrix technique which leads to a different, but equivalent, expression. This 
calculation is presented in appendix A, and it is expressed as 

'^(u h) 

' ~ trace[Tceii{(^,h)] = 2cos{dq), (30) 



Jl J2... Jn 



where Tceii{uj, h) is given by eq.(A.15). Although it is not shown in Appendix 

A, we have verified that '^{uj, h) depends on the square of the exchange con- 
stants, {Jf, J|, J^}. This means that the effect of the change of the signs 
of the J's is to introduce a shift of n/d in the spectrum wave-vectors, as can 
be seen in eq.(30). 

As it is also shown in Appendix A, for zero external field, there is no energy gap 
between the ground state and the first excited state, for n odd and arbitrary 
J's. However, for n even, the gap in the spectrum is only suppressed when the 
J's satisfy the condition 

J\J^...Jn-l = J2JA---Jn- (31) 

This latter result can also be obtained from the criticality condition for ID 
random Ising model in a transverse field [25] and the equivalence between 
XY chain and two decoupled transverse field Ising chains[26,27]. 

The excitation spectrum for a chain with n = 8, equal and different magnetic 
moments, is shown in Fig. 2. As can be seen, when the J's satisfy eq.(31) 
(continuous line) there is no gap at zero field, whereas a gap is present when 
eq.(31) is not satified (dot-dashed line). Although these results are presented 
for equal /x's, they are still valid when the /x's are different. 

The dotted line represents the spectrum for nonzero field, different ji's and 
the J's also satisfying eq.(31). As can be seen, the spectrum shifts and opens 
a gap at zero and boundary wave- vectors, which is a consequence of the non- 
comutativity, in this case, of the field term with the Hamiltonian. 
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3 The induced magnetization and isothermal susceptibility Xt 



Prom eqs.(3) and (22) we can express the local induced magnetization M^^, 



z 

m 



as 



Ml,m = /^m (Slm) ^ K,km'^1'kmnq,k - (32) 

q,k 

where the occupation number Ug^k is given by 

^^.^ = rr^' ^^^^ 

and the calculation can be done by considering H — H~ [20,21,22], since we 
are interested in the thermodynamic limit. 



We can define an average cell magnetization operator in the z direction, rf , 
as 

1 " 1 / 1\ 

= - E l^rnS[,m = " E ( 4,mCl,m " ^ ) , (34) 
m=l "' m=l ^ 

which corresponds to a generalization of the cell spin operator defined in the 
study of the alternating supcrlattice[18]. By using eqs. (32) and (34) we can 
write the induced magnetization per site as 

1 " / 1 i\ 

= i^l) = - E /^"^ hv E K,kmUq,kmng,k -7:], (35) 
^ m=l q,k ^ 



which can be written in the form 



E l^rnU*g^km'^q,km tauh ( 



f3£q,k 



= 22 ^imUlkm'^q,km tanh { ) . (36) 

The isothermal susceptibility can be obtained from the expression 



Xt = z^T- = \ E f^m—^Ug,kmtBnh - 

\^q,k,m \ 



n dh 2nN 
+ E /^m<fe„,^^^tanh (^y^J + 



+ 2 E ^^mUg,kmUq,kmSech I — - I —— \ . (37) 
q,k,m \ / ) 

At T — 0, where the system presents quantum transitions, the induced mag- 
netization, obtained from eq. (36) in the limit T — > 0, is given by 

= -^-^ E l^mUlkm1^g,kmSign {Sg^k) , (38) 
q,k,m 
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and from eq.(37) wc obtain xf?, which diverges at the critical fields he- We 
identify the largest critical field as hg, since for h ^ hg the induced magneti- 
zation is saturated. 

For identical magnetic moments, the average cell magnetization operator, rf , 
is proportional to the average cell spin operator, ^ Z]^=i Sf^^, and eqs (36) 
and (38) can be written, respectively, as 

= E .tanh (^) , (39) 

q,K,m \ / 

= E/^sign(%fc)' (40) 

q,k,m 

where /i^ — /i for any m. 

The results for and xf?, at T = 0, are presented in Fig. 3 as functions of 
the field /i, for a chain with n — 8 and identical /I's. The continuous line and 
the dashed line correspond to the cases where the exchange constants satisfy 
and do not satisfy the condition shown in eq.(31), respectively. As we have 
shown for the alternating superlattice [18], the magnetization also presents 
plateaus which are limited by critical fields he, where the isothermal suscep- 
tibility diverges, which correspond to quantum phase transitions induced by 
the field. The regions of plateaus, which we associate with disordered regions 
[18], correspond to the gaps in the excitation spectrum, and the critical fields 
are associated with the zero-energy mode with the wave-vectors q = and 
n/d. 

As it can also be seen, when there is no gap in the excitation spectrum at 
zero field (continuous fine), the zero magnetization plateau, which is present 
when there is a gap (dashed line), is suppressed, and consequently the total 
number of transitions induced by the field is n — 1. It should be noted that 
the local magnetization also presents plateaus and non-analytic behaviour at 
the critical fields. 

For different magnetic moments, the critical fields are obtained as in the pre- 
vious case, and the results for and Xt^ ^^^^ at T = 0, are shown in Fig. 4. 
The main difference from the previous case is the fact that the magnetization 
does not present plateaus, and the susceptibility Xt different from zero at 
both sides of the transition, although it diverges at one side only. The local 
magnetization presents similar behaviour as the total magnetization, and there 
is no suppression of the transition at zero field when there is no gap in the 
excitation spectrum (continuous line). The regions between two critical fields 
where the susceptibility Xt is finite correspond to the gaps in the spectrum. 
As in the case of equal magnetic moments, we associated these regions with 
disordered regions. 
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It should also be noted that for T > 0, all these transitions are suppressed by 
the thermal fluctuations. 



The exact expressions for the critical fields for n = 2, 3 and 4, when we have 
different magnetic moments, can be obtained from the analytic expressions for 
the spectra given in eqs.(23),(24) and (27), by considering Sq — and q — 
or Tr/d. The explicit results are 



cl 



IJ1 + -/2I 



c2 



\Jj_^_J2\ 

2VA*iAf2 



■, forn = 2 [14], 



(41) 



'9' + 2n 



cos 



hc3 — 2\JR'i 
where 



cos 



'9' + -k' 



, for n = 3, 



(42) 



9' — arccos 



Rf 

_ J^lls + J|/ii + J3V2 



2 ' with < ^73 < 7r/2, 
J1J2J3 



12/^1/^2/^3 



8/^1/^2/^3 ' 



and 



/iei = /ii+-\ K2 = h';++\ h,s = hi-+\ /ic4 = /ii" forn = 4, 
where 



(43) 



(44) 



/fl/t2/t3/t4 



±- 



1 



"2/fi/t2/f3/f4 
-4/ii/X2/i3/^4 



|^(/il/t2-''| + /il/^4-^2 + /i2/i3</4 + /^3/^4-^l)^- 



1/2 >| 1/2 



(45) 



For arbitrary n, the critical fields can also be obtained from the solution of 
the equation 
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For homogeneous magnetic moments we can also obtain from eq. 
T = 0, analytical expressions for the induced magnetization . For 
values of J's, the critical fields satisfy the relation hd > hd+i and the 
magnetization, for n = 2, is given explicitly by 



(38), at 
positive 
induced 



h>hci = hs, 
^ arccos A2 hc2 < h < hd, 
h < hc2, 



where 



A, 



2/^2/^2 _ 1 ( J2 + J, 



J1J2 



with hci and hc2 given in eq.(41), and for n = 3 by 



2 

f - ^ arccos(-A3) 
a 

6 

-^arccos A3 



h>hci = hs, 
hc2 < h < hci, 
hc3<h< hc2, 
h < hc3, 



where 



A, 



J1J2J3 

with hci, hc2 and hc3 given in the eq.(42), and finally for n = 4 by. 



2 

f - ^arccos A4 



4 ' 



h>hci = hs, 
hc2<h< hci, 
hc3<h< hc2, 



f - ^ arccos(-A4) hc4 < h < hcs. 



h < hci, 

where hd, hc2, h^ and h^ are given in eq.(44), and 



. 8/.^/.^ - 2 ( J| + J| + J| + Jl) ,,'h' {JfJl + JlJl) 

»-4 — ^ ^ ^ ^ h 



J1J2J3JA 



2J1J2J3J4: 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



For arbitrary n, between two critical fields, namely, hcj and hcj+i, with hd = 
hs, a general expression can be obtained for M? as a function of h and is given 

by 



Mf 



P3 



arccos 



nn 



n+ 



7+3 



A. 



(53) 
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where M^j — ii[n — 2{j — 1)] /2n is the magnetization in the upper plateau 
and An is given by 

A„ = ^^^. (54) 

2 Ji J2... Jn 



Since the model has azimuthal symmetry, for homogeneous magnetic mo- 
ments, the values of the induced magnetization per site at the plateaus, Mp, 
satisfy the quantization condition [28] 

n{iJ,/2 — Mp) = ij, X integer, (55) 

which is verified by the results shown above in eqs. (47-51). For different mag- 
netic moments, although the azimuthal symmetry is preserved, this condition 
is no more satisfied since the total magnetization operator does not commute 
with the Hamiltonian. However, it must be noted that eq.(55) is always sat- 
isfied provided we replace the induced magnetization per site by the average 
spin component in the field per site and assume = 1. 

For different /I's, when there are no magnetization plateaus, we can define an 
order parameter associated with each second order quantum transition given 
similarly to the case where we have identical /I's [17] as 

^\M^ - M^\, (56) 



where is the magnetization at the critical point, and from this we can 
show that the critical exponent (3, obtained from the scaling relation 

M'r^\h-hcf, (57) 



is equal to 1/2, and the exponent 7 associated with the isothermal susceptibility, 
obtained from 

x¥-\h-K\~\ (58) 

is also equals to 1/2. 

For the special case when the Jfs satisfy the special condition, eq( 31), the 
critical exponents j3 and 7 associated with the transition that occurs at h = 
are given by 1 and respectively. Therefore, the quantum transitions of 
the inhomogeneous model, apart from this special point, belong to the same 
universality class as the ones in the homogeneous model [29] . 
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4 Static and dynamic correlations (riTf^j.'j 

As it has been shown [20,21,22], in the thermodynamic hmit, the correlation 
function ^>S'f^m(t)5'f_j_j.^,(0)^ can be obtained from the expression 

_ Tr[eM-PH-)eMiH-t)Sl^eM-iH-t)Sf^r,n 

,^l,m{t)^l+r,m'{^)) " j.^ [exp(-/3i/-)] 

(59) 

and consequently we can express this correlation in terms of fermion operators. 
Therefore, the dynamic correlation between the effective spins in the field 
direction 

1 n 

(rf Wr,V(0)) = - E l^mf^m'{SUt)SUr,m'iO)), (60) 
^ m,m'=l 

can be written as 



T[{t)Ti\M) = Yl ^rn^ra' (cJ,m(^)Q,m(t)c^^,(0)Q+,,^.(0)) - 



m,m'=l 



(^)Q,m(^)) + (c;V,m'(0)Ci+r,m'(0))) + 



+ (61) 



2 

4 

and by using Wick's theorem [30], we obtain 
1 

{rnt)ri\M) = -I E l^ml^m' {cl^{t)ci+r,m'{^)) (Q,.n(t)qV,^.(0)) + (rf )^ . 

ft I -I 

m,m — 1 

(62) 

From eqs.(22), we can write the contractions 

(4,m(^)Q',m'(0)) = ^ 5^e^^<^('-'')ti;fc^ti,,fc^.n,,fcexp(i£,,fei), 

q,k 

[ci,m{M',mm) = ^ Ee"'''^'-''^^.,fem<,fen.' (1 " n,,k) exp(-i£,,,i) , (63) 

q,k 

and from this equation, we obtain 



Tf(t)Tf+^(0)) = -^ ^ IImf^rn'Ulkm1^g,km'Ug>,k'mU*g,^k,^, X 

X (1 - n,,,fcO e-^''^(«-^')e<^-'=-V,.')* + (rf , (64) 
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which in the hmit T" — > can be written as 



f V ly I 111 I 

q,k q',k m,m' 

{for £q,fc<0) {for V.fc'>0) 



The static correlation (Tir[_^j.J is immediately obtained from the dynamic one 
by making t — and we obtain, for finite T, 



'zz 
■Jl U+r 



1 



„2/\r2 ^ ^ r-inr''m'Uq^},^Uq^krn'Uqi ^k'mUqi ^^1^1 'X 
I L ly J I 111 
q,m,k q',m,k' 



xnq,kil-nq,,k')e-'''^'^-''^ + {Tn\ 



(66) 



and 



^z z 

'I '1+1 



n- 



E EE t^mi^m'Uqkm'^q,km'X 

q,k q\k' m,m' 

{for £,,fc<0) (for £,/,fc/>0) 



^''^q',k'mUqi^k'm'^ 



-idr{q—q') 



+ {riT , 



(67) 



in the limit T — > 0. 

For different //'s, we present in Fig. 5, for a chain with n = 8, at T = 0, the 
static correlation as a function of the distance between cells for field values 
above, below and at a critical point. As expected, for fields below the critical 
one, which in this case corresponds to the region which contains the diver- 
gence of the isothermal susceptibility, the correlation presents an oscillatory 
behaviour. As in the alternating superlattice [17], the period of the oscilla- 
tion corresponds to the correlation length and goes to infinity as we approach 
the critical field, and there is no oscillation in the correlation for fields in the 
disordered region. 

The real and imaginary parts of the dynamic correlation as function of the 
time, for a chain with the same set of parameters of Fig. 5, are presented in 
Fig. 6. Differently from the case where we have identical fi's, in the region 
between two critical fields, where the isothermal susceptibility is finite, the 
dynamic correlation varies with the field, and for h ^ hs is time-independent 
and equal to (rf )^. 
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5 Dynamic susceptibility x^^i'^) 

Introducing the time Fourier transform in tlic dynamic correlation (^Tf (t)T;+r(0)^, 

(^^^^+^). ^ ^ /I (^''(^)^'V(O)) e^^'rft, (68) 
we obtain from eq. (64) the result 

ltl\ J 111 I 

X (1 - n,,,k') e-"'<'^-^'^5 (a; - e,,, - e,. ,0 + (rf 5(a;). (69) 
By considering the spacial Fourier transform in the field direction, 

(-^-.l-E(-^f..>.e--, (70) 

r 

we obtain the correlation as function of uj and g, in the form 



T„ T 



1 



q' k,mk',m' 

X ''^g'.A: (1 — nqi_q^ki) 5 {uJ + Eqi^k — ^q'-q,k') + 



{t^Y 6{uj)5q,,. (71) 

From this expression we can determine the dynamic susceptibility Xq^{'^)-i 
which is given by [31] 



Xr((^)^-27r((r-r!,)), (72) 
where the Green function Ut^]t^X) is obtained from 



u — u' 



duo', (73) 



and we obtain 



Xr('^)-~:^^XX X] l^mfJ'm'U*g>^km'^q',km'Uq'-q,k'mU*q'-q,k'm' >< 

q' k,mk',m' 

^ nq',k - nq^-q,k' ^^^^ 
^ + £q',k — ^q'-q,k' 
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where we have used the identity 

1 - nq,k = ef^'^'^Ug^k- (75) 

We can also write xTi^) the form 



q' k,mk,m' 

^ tanh(/3£g/^fc/2) - tanh(/?£g/_g^fc72) ^^^^ 

<^ + £q',k — £q'-q,k' 

and from this obtain the static susceptibihty Xo^(O) which is exphcitly given 

by 

Xo^(O) = 22/V ^ ^ ^ ^m^m'Ul, ^km^qi ,km'Uq' ,k'mU*qi ^k'rn'^ 
q' k,k' m,m' 
kj^k' 

tanh(/3£g/,fc/2) - tanh(/3£g/^fc//2) ^ 

£q',k — £q',k' 

+ ;^4]V^ ^ /^m/^m'<,Jfc,„V,fcm'V,fem<',Jfc,„'Sech2(/3£g/,fe/2). (77) 
q'fc m,m' 

As in the uniform model [32], it can be shown from the previous expression 
that Xo^(O) is equal to the isothermal susceptibility Xt given in eq.(37). 

The real and imaginary parts of x^^(a;), as usual, arc obtained by considering 
Xl^i'^ ~ i^) in the limit e ^ in eq.(76) and are given by 



R^Xq^ i'^) — 2ni\r^'^' X] l^rnl^rn'U*qt^km^q'Mrn'Uqi-q,k'mKi'-q,k'rn' ^ 
q' k,mk,m' 

tanh(/3£g/,fc/2) - tanh(/3£g/_g,fc72) ^^^^ 

<^ + £q',k — £q'-q,k' 

where P denotes Cauchy principal value, and 



Imxl\uj) = ^ ^ ^ l^mlJ'm'Ki',km^q',km'Uq'-q,k'rnU*qi-q,k'm' ^ 

q' k,m k',m' 

x(tanh(/3£q.,jfc/2) - tanh(/3£g/_g,jfc72)) x 

xS{u; + eq>^k-eq'-q,k')- (79) 

At T = 0, these results, as function of ur, are shown in Fig. 7, for n = 4 and 
different wave-vectors for a field in the disordered region. For any wave- vector 
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there are the same number of bands in the imaginary part of the susceptibihty, 
and this number depends on the size of the unit cell and on the values of the 
parameters. 

In Fig. 8 it is shown the static susceptibility, x^^(O), at T = 0, for the same 
lattice parameters of Fig. 7, as a function of the field. For q = 0, since it is 
identical to the isothermal susceptibility, it diverges at the critical fields. The 
singularities at non-zero and at the Brillouin zone boundary wave- vectors can 
be associated with critical points, whereas the ones for different wave-vectors 
are related to oscillations of the spin correlations and can also be associated 
with the unstable critical points present in the study of the system within the 
real space renormalization group [33]. 



6 T=0 isothermal and dynamic susceptibilities xW Xq^i'^) 

At T = 0, and for h> hs, we can obtain from the eq.(B.8) the time Fourier 
transform of dynamic correlation (Si^^{t)Sf_^_^^^,{0)\ , which is given by 



{SlA,m')^ = ^ r exp(^a;t) (SlMS^+rM^)) dt, (80) 



oo 



27r 

and can be written as 

^l,m^l+r,m')^ = ^ Z ^M-'iQdr)ulk^Ugkm'S{uJ + (81) 

q,k 

Following Gongalves [34], we can write the local dynamic susceptibility x^^(to', r, m, m') 
as 



X ioj, r,m,m') ^ - hm / ^— '-^du', (82) 

since < 0, and by using eq.(81) we obtain 

X u;,r,m,m = — . (83) 

The susceptibility associated with the cell effective spin in the x direction is 
given by 

^^^(a;,r) = X''''{'^,r,m,m'), (84) 

m,m' 

and by using eq.(83) we can write 

= _(zf: J2 (-l)'"- """'"'^*''"''"""'''""' (86) 

4Ar ^ / ^ UJ + Eak 

q,K,m,m' 'i''^ 
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and from this the final result 



u r 4 (a; + fc ) 



where 

q for n even, 



(87) 

0' — I for n odd. 

From eq.(86) by considering the x^^i^o — ie) in the hmit e — > we can obtain 
immediately the real and imaginary parts which are given by 



1 \m— m' * 

k,m,m V y 



where P denotes Cauchy principal value, and 

7mxfH = -^ E i-'^r'"''K',km^i'>'^m'S{u; + e,,,,). (89) 

We can also obtain, in the limit considered, the local isothermal susceptibility 
from the expression [35] 

(r, m, m') = lim (^l%(-^A)5,%,^,(0)) dX, (90) 

which, from eq.(B.8), can be written as 

Prom this result we can immediately conclude that the isothermal susceptibil- 
ity) XW — Sr,m,m' xWi''^: equal to the static one, Xq^{^): in the limit 
g ^ 0. 

For equal magnetic moments, n=8 and the exchange parameters of Fig. 3(a), 
the real and imaginary parts of the dynamic correlations (^Si ,^{t)Sf_^_^j^,{0)'j 

and ('S'f^(t)S'f_^^^,(0)^ are shown in Figs. 9 and 10 respectively, for various 
magnetic fields. For large t, these correlations factorize and tend to zero, since 
the average (5";%) ^ 0. 

Due to the simple form of the dynamic susceptibility, eq.(86 ), it presents poles 
of order one at the excitation energies. This means that the imaginary part 
corresponds to a set of delta functions for a given wave-vector. 

In Fig. 11 we present the static susceptibility Xg^(a;), at T = 0, as a function 
of the field for a chain with n — 8. For the ferromagnetic case shown in Fig. 
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11(a) it diverges at the critical field for q = only. On the other hand, if we 
change the sign of a single exchange interaction, it diverges for q at the zone 
boundary, namely, q = t:/8. This is shown in Fig. 11(b), and means that the 
system in terms of cells has an antiferromagnetic behaviour. Although it is 
not presented, it can be shown that the order ferro (antiferro) is associated to 
the sign + (-) of the product of the exchange interactions within the cells. 



7 Conclusions 

We have considered in this work the isotropic XY model on the inhomogeneous 
periodic chain with cells, n sites per cell. The exact solution has been 
obtained in the general case, where we have n different exchange constants 
and magnetic moments, and for arbitrary temperatures. 

At T = 0, for equal magnetic moments, we have shown that the induced 
magnetization, as a function of the field, presents plateaus which satisfy the 
quantization condition shown in cq. (55), and that the isothermal susceptibility 
Xt is zero within the plateaus and diverges at the limiting fields which char- 
acterize the multiple second order quantum transitions. On the other hand, f 
or different magnetic moments, although the system presents multiple quan- 
tum transitions, the induced magnetization plateaus are suppressed and the 
isothermal susceptibility, which also diverges at the critical points, is always 
different from zero. In both cases, the disordered phases correspond to the 
gaps in the excitation spectrum and to the regions where Xt finite, which 
are identical to the plateaus for the case of identical magnetic moments. For 
n = 1,2,3, explicit analytical expressions are presented for the critical fields, 
in the general case, and for the induced magnetization when we have identical 
magnetic moments. 

The number of transitions, for different fi's, is always equal to n, even for the 
special case where the exchange constants, for n even, satisfy the condition 
presented in eq.(31), where at zero field there is no gap between the ground 
state and the first excited state. It should be noted that in this case, Xt 
diverges at both sides of the zero field transition and that, for identical /I's, 
this transition is suppressed. 

The critical exponents are identical to those of the uniform model, a = 1/2, 
P — 1/2, 7 = 1/2, and for the special point at zero field a = 0, /3 = 1, 7 = 0, 
and naturally they satisfy the Rushbrook relation a+ 2/3+ 7 = 2 (see e.g. 
ref.[36]). It should be noted that in these transitions a = 7 [37]. 

The correlation length, which corresponds to the period of the oscillation of 
the static correlation, (rf r;!)_^\ , as in the homogeneous chain and alternat- 
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ing superlattice [17], diverges at the critical points irrespective of the val- 
ues of the magnetic moments. On the other hand, the dynamic correlation 
(t[ {t)T[^^{Q)'j , which is independent of the field in the plateaus regions for 
equal magnetic moments, is field dependent when we have different magnetic 
moments. 

The isothermal susceptibility Xt equal to the the static one, Xo^(0)i 
any temperature. Independently of the value of the field and of the magnetic 
moments, the imaginary part of the dynamic susceptibility Xq^i'-^) presents 
several bands whose number does not depend on q. As in the alternating 
superlattice, the discontinuities of the imaginary part at the band edges cor- 
respond to singularities in the real part. Associated with the critical behaviour 
are the descontinuities which occur at zero or at the Brillouin zone boundary 
wave-vectors, whereas those that occur for other values of the wave-vectors 
are related to unstable critical points [33] . 

At T = and for fields greater than the saturation field, we obtained the 
dynamic correlation <^'S'f m(t)5'f_,_^„j/(0)^ and the isothermal and wave- vector 
dynamic susceptibilities in the x direction. We have shown that the isother- 
mal susceptibihty Xt^ is identical to the static one Xo^(O), and that in terms 
of the magnetization of the cells the system presents a ferromagnetic or an- 
tiferromagnetic order depending on the sign of the product of the exchange 
interactions within the cells. 
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A Transfer matrix technique: excitation spectrum 

The excitation spectrum of Hq, 



m=l 

n-1 



J n 



71=1 ^ 



A ' A 

[^S,n^?,iexp(-igrf) + Aj^iAg,„exp(igd) 



(A.l) 



given in eq.(19), can also be obtained by using a transfer matrix method. By 
introducing the equation of motion for the operators Ag^^n 



(A.2) 



and assuming the time-evolution given by 



Aq,m{t) = Aq,^exp(-ia;t), 



(A.3) 



we obtain the system of equations 



. , Ji . Jncxp{iqd) 

[U + = -y^<z,2 ^g,n, 



/ I U\ A Jm—l 4 Jm—2 a 

[LJ -+- IJ'm-l'^)^q,m-l — ^ ^q,m ^ ■^q,m-2, 



(a; + IJ,n_ih)Aq^n-l — '^^q,n ^^(j,n-2j 

. . , X . Jnexp{-iqd) J^-i . 

(a; + UnhjAq^n = ^ A^^i —Ag^^^i, 
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1.5 




Fig. 2. Excitation spectrum for n = 8, Ji = 1, J2 = 2, J3 = 1.5, J4 = 5/3, J5 = 2/3, 
Je = 3/5, J7 = 2, Jg = 1 for /xi = .... = /xg = 1 and /i = (continuous line) and for 
fj,i = 3, fi2 = 4, Us = = ^17 = 1-5, /i4 = 2.5, HQ = fis = ^ and /i = 0.2 (dotted 
line), and for Ji = 2, J2 = 1, J3 = 1.5, J4 = 5/3, J5 = 2/3, Jg = 3/5, J7 = 2, 
Jg = 1, Hi = .... = yUg = 1 and h = (dot-dashed line). 



23 



0.50 



(a) 



-Hr 



0.25 



/ 



0.00 



/ 



0.0 0.2 0.4 0.6 0.8 

h 



1.0 1.351.401.45 



1.5 



^1 



zz 



1.0 



0.5 



0.0 



(b) 



u 



^ /I 



0.0 0.2 0.4 0.6 0.8 1.0 1.351.401.45 



h 



Fig. 3. (a) Magnetization and (b) isothermal susceptibility Xt^^^ T = and 
for n = 8, /ii = .... = = 1; as functions of the uniform field for Ji = 1, J2 = 2, 
J3 = 1.5, J4 = 5/3, J5 = 2/3, Je = 3/5, J7 = 2, Jg = 1 (continuous line) and 
= 2, J2 = 1, J3 = 1.5, J4 = 5/3, J5 = 2/3, Je = 3/5, J7 = 2, Jg = 1 (dashed 
line) . 
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(a) 




0.0 0.2 0.4 0.6 0.90 0.91 

h 



Fig. 4. (a) Magnetization and (b) isothermal susceptibility Xt^ at T = and 
n = 8, jUi = 3, jLt2 = 4, = /Lt5 = /Lt7 = 1.5, jLt4 = 2.5, jLte = /xg = 1 , as functions 
of the uniform field for Ji = 1, J2 = 2, J3 = 1.5, J4 = 5/3, J5 = 2/3, Jq = 3/5, 
J7 = 2, Jg = 1 (continuous line) and Ji = 2, J2 = 1, J3 = 1.5, J4 = 5/3, J5 = 2/3, 
Je = 3/5, J7 = 2, Jg = 1 (dashed line). 
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Fig. 5. Static correlation function < rfrf^^^ > as a function of r (distance between 
cells), at T = 0, for n = 8, Ji = 1, J2 = 2, J3 = 1.5, J4 = 5/3, J5 = 2/3, Jq = 3/5, 
J7 = 2, Jg = 1, fj,i = 3, /i2 = 4, fj,^ = = fj,y = 1.5, //4 = 2.5, /ig = = 1 for 
values of the field near and at the critical field. 
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Fig. 6. The real and imaginary parts of the correlation function ("rf (t)Tj^i(O)) , as 
functions of time, at T = 0, for n = 8, Ji = 1, J2 = 2, J3 = 1.5, J4 = 5/3, 
J5 = 2/3, Je = 3/5, J7 = 2, Js = 1, m = 3, /X2 = 4, = = ij,^ = 1.5, /X4 = 2.5, 
= A^s = 1, for values of the field near and at the critical field. 
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Fig. 7. The real and imaginary parts of the dynamic susceptibility in the field 
direction, Xq^{^)^ at T = as a function of frequency for n = 4, Ji = 1, J2 = 2, 
J3 = 1.5, J4 = 5/3, /xi = //3 = 1, = A*4 = 0.5 and /i = l,for different values of q. 
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0.4 




Fig. 8. Static susceptibility in the field direction, Xq^i^)^ at T = 0, as a function of 
the field for n = 4, Ji = 1, J2 = 2, J3 = 1.5, J4 = 5/3, fj,i = 1, fX2 = 0.5, /X3 = 1, 
H4 = 0.5 and different values of q. 
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Fig. 9. The real and imaginary parts of the correlation function (^Sf^{t)Sf^{0)j , as 
functions of time, at T = 0, for n = 8, Ji = 1, J2 = 2, J3 = 1.5, J4 = 5/3, J5 = 2/3, 
Jq = 3/5, J7 = 2, Js = 1, m = .... = /X8 = 1, for values of the field at and at above 
the saturation field hs. 
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Fig. 10. The real and imaginary parts of the correlation function (^Sij^{t)Si^j.-^{0)j , 
as functions of time, at T = 0, for n = 8, Ji = 1, J2 = 2, J3 = 1.5, J4 = 5/3, 
J5 = 2/3, Jq = 3/5, J7 = 2, Jg = 1, /xi = .... = yug = 1) at the saturation field hg. 
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(b) 

3000 n , 




Fig. 11. Static susceptibility Xq^i.^) ^ ^ function of the field, at T = 0, for n = 8, 
Hi = 3, /i2 = 4, = /X5 = /xy = 1.5, /X4 = 2.5, jLte = /xg = 1 , and different values 

of q. (a) Corresponds to Ji = 1, J2 = 2, J3 = 1.5, J4 = 5/3, J5 = 2/3, Jc, = 3/5, 
J7 = 2, Jg = 1 and (b) Ji = -1, J2 = 2, J3 = 1.5, J4 = 5/3, J5 = 2/3, Je = 3/5, 
J7 = 2, Jg = 1. 
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which is equivalent to the eigenvalue equation shown in eq.(14). Prom this set 
of equations we can write the matrix equation 



(a 

A 



T^(a;,/i) 



'■q,m 



, for m = 2, 3, n — 1, 



where 



Tm{uj, h) 



_n / ijJ+JMnh\ Jm-1 

Jm Jrr 



and 



(A.5) 



(A.6) 



A, 



= Ti(a;,/i) 



A ,\ ( 



A, 



I Aq n-l , 



for m — 1, 

, for m = n 



with T 1 and T „ given by 



Ti{u;,h): 



2(^) -^exp(^dg) 



V 



V 



T„(a;, h) = exp{idq) 



y 



^exp{—idq) y 
which satisfy the result, 

Ti(ci;, h)Tn{w, h) = exp(mg)Ti(cc;, /i)T„(a;, h), 
with Ti and T„ given by 



T„(a;, h) = 



Jl I Jl 



Of UJ+IJ,„h \ Jn-l \ 

Jn > J, 



(A.7) 

(A.8) 



(A.9) 
(A. 10) 



(A.ll) 



(A.12) 
(A. 13) 
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Prom eqs.(A.5-A.13) we can write 



Tcell{(jJ,h) 



' \=exp{-iqd)\ ' I = exp(-ig(i)*5,i,2, (A.14) 



where 



T^eiiii^, h) = Ti(a;, /i)T„(a;, /i)T„_i(a;, /i)T„_2(a;, h) . . . ¥3(0;, h)T2{u, h). 

(A.15) 

Eq.(A.14) shows exphcitly that exp(— ign) is an eigenvalue oiTceiii^^^h) cor- 
responding to the eigenvector ^E'g,i,2- Since Tcen(t^j h) does not depend on g,we 
can obtain immediately the second eigenvector of Tceii{oJ-, h), and the respec- 
tive eigenvalue, from this equation by introducing the transformation q —>■ —q, 
which gives^_5^i_2 and exp(— igo?) respectively. These results are consistent 
with the fact that det[Tceii{i^, h)] — 1, and we can write finally the equation 

trace[Tceii{uj, h)] — 2cos{dq), (A. 16) 

whose solution will give the excitation spectrum. 

In the absence of the external field, for n odd, we can show that the previous 
equation can be written in the form 



f{uj)xu;^2cos{dq), (A. 17) 

where f{uj) is a polynomial function of degree n — 1. This means that q = 7r/2d 
is a zero-energy mode irrespective of the values of J's. 

On the other hand, for n even and zero external field, eq.(A.16) can be written 
as 



where g{u!) is a polynomial of degree n — 2. From this result we can con- 
clude that g = is a zero-energy mode provided the condition JiJ^...Jn-i = 
J2JA---Jn is satisfied. Therefore, under this condition, there is no energy gap 
between the ground state and the first excited state in the absence of external 
field. 



JoJa-'-Jti. 



+ 



Jl J3... Jn-1 



JlJs...Jn-l J2Ja---Ji 



2cos(dg), (A. 18) 
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B Dynamic correlation {Sf^Jt)Sf^r,m'{^)) 



The dynamic correlation function \^Sf^{t)Sf_^^„^,{0)j , in terms of the hamil- 
tonians H^, in the thermodynamic hmit, can be written in the form [20,21,22] 



Tr 



exp{-(3H-) exp{zH-t)Sl^ exp(-^if+t)5f+. 



Tr [exp{-l3H-)] 



(B.l) 

which is very difficult to calculate, since H and do not commute. How- 
ever, at T = 0, for external field greater than the saturation field {h ^ hs), 
and H~ have identical ground state, namely, 

^J{®\nim) (B.2) 

l,m 

where nim = 1, V m, and 

H~ |$o> = H+ |$o) = |$o) ■ (B.3) 
Then, in this limit, the dynamic correlation ('S'^,TO(^)'^i+r,m'(0)) given by 
{Sl^it)SUr,A^)) = {^o\exp{tHn)Sl^exp{-tH-t)St^^^^, \%) . (B.4) 

Prom the eqs.(2) and (3) we can express Sf,^ in terms of fermion operators in 
the form 



m'-l 



i=l k=l k'=l 



r,k' 



and substituting this result in eq.(B.4) we obtain 



(B.5) 



l,m 



which can be written as 



(-1) 



nr+m'—l 



exp{iEot){-ir-' ($0| {elm + Cl,m) X 
X exp{-iH-t) [c\^r,m' + Cl+r,m') |$o) , (B.6) 



\nr+Tn' —m 



{^o\exp{iH i)[cl^„ + ci,jexp(-ii7 t) 



(B. 
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Finally, by using eq.(63), we obtain 

q,k 

(B.8) 

which reduces to the known result for the homogeneous chain [38]. 

Prom the previous equation we can conclude immediately that the static cor- 
relation {SlrnSf+r,m') ^qual to 5g,o/4. 
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